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Abstract 

The cohesive energy of condensed matter in strong magnetic fields is a fundamen- 
tal quantity characterizing magnetized neutron star surfaces. The cohesive energy 
refers to the energy required to pull an atom out of the bulk condensed matter at 
zero pressure. Theoretical models of pulsar and magnetar magnetospheres depend 
on the cohesive properties of the surface matter in strong magnetic fields. For exam- 
ple, depending on the cohesive energy of the surface matter, an acceleration zone 
("polar gap") above the polar cap of a pulsar may or may not form. Also, con- 
densation of the neutron star surface, if it occurs, can significantly affect thermal 
emission from isolated neutron stars. We describe our calculations of the cohesive 
property of matter in strong magnetic fields, and discuss the implications of our 
results to the recent observations of neutron star surface emission as well as to the 
detection/non-detection of radio emission from magnetars. 

Key words: stars: pulsars, stars: neutron, stars: magnetic fields, radiation 
mechanisms: non-thermal 



1 Introduction 



Recent observations of neutron stars have provided a wealth of information 
on these objects, but they have also raised many new questions. For example, 
with the advent of x-ray telescopes such as Chandra and XMM-Newton, de- 
tailed observations of the thermal radiation from the neutron star surface have 
become possible. These observations show that some nearby isolated neutron 
stars (e. g., RX J1856.5-3754) appear to have featureless, nearly blackbody 
spectra ( Burwitz et al.l 20031 ; van Kerkwijk fc Kaplan! [20061 ) . Radiation from 
a bare condensed surface (where the overlying atmosphere has negligible opti- 
cal d epth) has been invoked to explain this nearly perfect blackbody emission 



e.g.. lBurwitz et al.ll2003l : Mori &: Rudermanll2003l : Ivan Adelsberg et al.ll2005 
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Turolla et al.ll2004l ; iPerez-Azorin et al.ll2006l ). However, whether surface con- 
densation actually occurs depends on the cohesive properties of the surface 
matter. 



Equally puzzling are the observations of the recently-discover ed anomalous x- 
ray p ulsars (AXPs) and soft gamma-ray repeaters (SGRs) 
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Woods fc Thompson 



20051 for a review). Though these stars are believed to be magnetars, neu- 
tron stars with extremely strong magne tic fields (B > 10 14 G), they mostly 
show no pulsed radio emission (but see Camilo et al. 20061 ) and their x-ray 
radiation is too strong to be powered by rotational energy loss. By con- 
trast, several high-B radio pulsars with inferred surfac e field strengths sim- 



ilar t o those of magnetars have been discovered (e.g., iKaspi fc McLaughlin 



20051 ). A deeper understanding of the distinction between pulsars and magne- 



tars requires further investigation of the mechanisms by which pulsars and 
magnetars radiate and of their magnetospheres where this emission origi- 
nates. Theoretical models of pulsar and magnetar magnetospheres depend on 
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2006). For example, depending on how strongly bound the surface matter is, 



an acceleration zone ("polar gap") above the polar cap of a pulsar may or 
may not form, and this will affect pulsar radio emission and other high-energy 
emission processes. 



The cohesive property of the neutron star surface matter plays a key role in 
these and many other neutron star processes and observed phenomena. The 
cohesive energy refers to the energy required to pull an atom out of the bulk 
condensed matter at zero pressure. For magnetized neutron star surfaces this 
cohesive energy can be many times the corresponding terrestrial value, due to 
the strong magnetic fields threading the matter. 



In two recent papers (IMedin fc Lai l2006al Jbl. hereafter ML06a,b), we calculated 
the cohesive energy for H, He, C, and Fe surfaces at field strengths between 
B = 10 12 G to 2 x 10 15 G. We now wish to investigate some implications 
of these calculations. This paper is organized as follows. In §11 we briefly 
summarize the cohesive energy results of ML06 used in this paper. In §111 we 
discuss the possible formation of a bare neutron star surface. The conditions for 
the formation of a polar vacuum gap in pulsars and magnetars are presented in 
§IV. We summarize our results in §V. A more detailed/comprehensive study 
of these and related issues will be presented in a future paper (Medin & Lai 
2007, in preparation). 
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2 Numerical Calculations of Condensed Matter in Strong Mag- 
netic Fields 



It is well-known that the properties of matter can be drastically modified by 
strong magnetic fields. The natural atomic unit for the magnetic field strength, 
Bo, is set by equating the electron cyclotron energy fiwBe = h(eB/m e c) = 
11.577 B12 keV, where B\i = B/(10 12 G), to the characteristic atomic energy 
e 2 /a = 2 x 13.6 eV (where a is the Bohr radius): 

£> = I^ = 2.3505 xl0 9 G. (1) 



For b — B/Bq ^ 1, the usual perturbative treatment of the magnetic effects 
on matter (e.g., Zeeman splitting of atomic energy levels) does not apply. 
Instead, the Coulomb forces act as a perturbation to the magnetic forces, and 
the electrons in an atom settle into the ground Landau level. Because of the 
extreme confinement of the electrons in the transverse direction (perpendicular 
to the field), the Coulomb force becomes much more effective in binding the 
electrons along the magnetic field direction. The atom attains a cylindrical 
structure. Moreover, it is possible for these elongated atoms to form molecular 
chains by covalent bonding along the field direction. Interactions between the 
linear chains can then lead to the formation of three-dimensional condensed 



matter flRudermanlll974j ; iRuder et al.lll994l ; lLaill200lh . 



Previous calculations of the c ohesive energy of on e- di mensional chains and con - 
densed matter were done by IJoned (119851 . Il986l ) and iNeuhauser et al.l (119871 ) . 
These earlier calculations adopted some crude approximations (e.g., band 
structure was not properly treated by Neuhauser et al.) and gave somewhat 
conflicting results. These calculations were also restricted to moderate neutron 
star field strengths {B = 10 12 -10 13 G). 

Our new calculations (ML06a,b) are based on density functional theory. ML06a 
summarizes our calculations for various atoms and molecules in magnetic fields 
ranging from 10 12 G to 2 x 10 15 G for H, He, C, and Fe, representative of the 
most likely neutron star surface compositions. Numerical results of the ground- 
state energies are given for Hat (up to N = 10), He^v (up to N = 8), Cn (up 
to iV = 5), and Fejv (up to iV = 3), as well as for various ionized atoms. 
ML06b summarizes our calculations for infinite chains for H, He, C, and Fe in 
that same magnetic field range. For relatively low field strengths, chain-chain 
interactions play an important role in the cohesion of three-dimensional (3D) 
condensed matter; an approximate calculation of 3D condensed matter is also 
presented in ML06b. Numerical results of the ground-state and cohesive en- 
ergies, as well as the electron work function and the zero-pressure condensed 
matter density, are given for H m and H(3D), He^ and He(3D), C^ and C(3D), 
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Fig. 1. Cohesive energy Q s and molecular energy difference AEn for iron (where 
AEn = En/N — Ei) as a function of the magnetic field strength. The symbol Qoo 
represents the cohesive energy of one-dimensional chains, and I± and I2 are the first 
and second ionization energies of the Fe atom. 

and Feoo and Fe(3D). 

Some numerical results from ML06a,b are provided in graphical form in Figs. [1] 
and [2j Figure [1] shows the cohesive energies and molecular energy differences 
A£jv = En/N — Ei for Fe, where E\ is the atomic ground-state energy and 
E N is the ground-state energy of the Fe^ molecule. Fig. [2] shows the electron 
work functions for He, C, and Fe, as a function of field strength. 



3 Condensation of Neutron Star Surfaces in Strong Magnetic Fields 



As seen from Fig. [H the cohesive energies of condensed matter increase with 
magnetic field. We therefore expect that for sufficiently strong magnetic fields, 
there exists a critical temperature T crit below which a first-order phase tran- 
sition occurs between the condensate and the gaseous vapor. This has been 
investigated in detail for hydrogen surfaces (see lLai fe Salpeterlll997l ). but not 
for other surface compositions. Here we discuss the critical temperature values 
for the phase transitions of Fe surfaces. 
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Fig. 2. The work function as a function of B, for H, He, C, and Fe infinite chains. 

A precise calculation of T cr i t is difficult to obtain at present. We can get an 
estimate by considering the equilibrium between the condensed phase (labeled 
"s" ) and the gaseous phase (labeled "g" ) in the ultrahigh field regime (where 
phase separation exists). The gaseous phase consists of a mixture of free elec- 
trons and bound ions, atoms, and molecules. For a given temperature, phase 
equilibrium occurs when the total pressure of the gas equals that of the con- 
densate, and when each species in the gas is in chemical equilibrium with the 
condensate. To simplify the calculation we have assumed that the vapor pres- 
sure is sufficiently small so that the deviation from the zero-pressure state of 
the condensate is small; this is justified when the saturation vapor pressure 
P sat is much less than the critical pressure P crit for phase separation, or when 
T < T crit . 

We have calculated the equilibrium densities of the various species in the 
gaseous vapor, as a function of temperature. Some results are shown in Fig. [31 
for magnetic field strengths of B = 10 13 and 5 x 10 14 G. The critical tem- 
perature T crit , below which phase separation between the condensate and the 
gaseous vapor occurs, is determined by the condition p s ~ p g . Using the values 
for E N /N (the energy per atom in the Fe^ molecule), E s (the ground-state 
energy of the condensed Fe), and E n+ (the energy of the Fe n+ ion) presented 
in ML06, we find T crit ~ 6 x 10 5 , 7 x 10 5 , 3 x 10 6 , 10 7 , and 2 x 10 7 K for 
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Fig. 3. The densities of various species and the total gas density in the saturated 
vapor of a condensed Fe surface, as a function of temperature, for B\2 = 10 (left 
panel) and 500 (right panel) . The horizontal lines give the condensed matter density 
p s . The curves labeld by a p g x Z; n t" give the vapor density when the interal partition 
functions of the atoms and molecules are included, while the other curves do not 
include the internal partition functions. The critical temperature is set by p g p s . 

B\2 = 5, 10, 100, 500, and 1000. In terms of the cohesive energy, we have 

kT CTit ~0.08Q S . (2) 



Note that although our results for the vapor densities were derived using 
P g "C p s , while the condition for phase separation requires p g ~ p s , we may 
still use the results to obtain an estimate of T cr i t . The vapor density drops very 
rapidly with decreasing temperature, such that when the temperature is below 
?crit/ 2 (for example) this density is much less than the condensation density 
and phase transition is unavoidable. When the temperature drops below a 
fraction of T cr i t , the vapor density becomes so low that the optical depth of 
the vapor is negligible and the outermost layer of the neutron star then consists 
of condensed iron. The radiative properties of such condensed phase surfaces 
have been studied using a simplified treatment of the condensed matter (see 



van Adelsberg et al.l 120051 and references therein) . 



4 Vacuum Gap Accelerators of Pulsars and Magnetars 



A rapidly rotating, strongly magnetized neutron star induces a strong electric 
field with a component along the magnetic field lines of the star. The plasma 
surrounding the neutron star is assumed to be an excellent conductor, such 
that the charged particles in the plasma move to screen this parallel component 
of the electric field. The particles fill the magnetosphere with a charge density 
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of approximately (IGoldreich fc Julian! Il969l ) 

n b 1 

PGJ = 



2vrc 1-ntrl/c 2 



(3) 



where O is the rotation rate of the neutron star. 



Charged particles traveling outward along the open field lines in the polar 
cap region of the neutron star magnetosphere will escape beyond the light 
cylinder. To maintain the required magnetosphere charge density these par- 
ticles will have to be replenished by the neutron star surface. If the surface 
temperature and surface cohesive strength of a neutron star are such that the 
required particles are tightly bound to the stellar surface, those regions of the 
polar cap through which the charged particles are escaping will not be replen- 
ished. A vacuum gap will form in those region s, starting at the stellar surface 



and moving outward at a velocity v ~ c (e.g.. iRuderman fe Sutherland! 11975 



Cheng fc Rudermanlll975l : lUsov fc Melrose! Il996l : iGil et al.M2003h . Obviously, 
the existence of this vacuum gap depends on the cohesive energy and electron 
work function of the condensed neutron star surface. 



We note that in the absence of a vacuum gap, a polar gap acceleration zone 

based on the space-charge-limited flow may still develop (e.g-. lArons fc Scharlemann 
19791 : lHarding fc Muslim^ Il998h . 



4-1 Thermal Emission/Evaporation of Charged Particles 



For neutron stars with O • B p > (B p is the magnetic field at the polar 
cap), electron emission is relevant. The emission rate for electrons from the 
condensed surface can be calculated by assuming that these electrons behave 
like a free electron gas in a metal. The energy barrier they must overcome 
is the work function of the metal. The equilibrium electron charge density 
(balance between electron loss to infinity and supply by the surface) above 
the polar cap can be written as 

Pcjexp (C e - 4>/kT) , 0>C e /cT, 
Pe ^ < (4) 
Pgj, < C e kT , 



where is the electron work function and the parameter C e ~ 30 depends 
on magnetic field strength and temperature. The electron work function was 
calculated in ML06b and is depicted in Fig. [2J 
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For neutron stars with f2 ■ B p < 0, ion emission is relevant. The emission rate 
for ions from the condensed surface is more complicated. Unlike the electrons, 
which form a relatively free-moving gas within the condensed matter, the ions 
are bound to their lattice sites. To escape from the surface, the ions must 
satisfy three conditions. First, they must be located on the surface of the 
lattice. Ions below the surface will encounter too much resistance in trying 
to move through another ion's cell. Second, they must have enough energy 
to escape as unbound ions. This binding energy that must be overcome will 
be labeled Eb- Third, they must be thermally activated. The energy in the 
lattice is mostly transferred by conduction, so the ions must wait until they are 
bumped by atoms below to gain enough energy to escape. The characteristic 
frequency with which that occurs is the lattice vibration frequency. 

From these three requirements for the ions, we can obtain the equilibrium ion 
charge density above the polar cap: 

( p G j exp (d - E B /kT) , E B > QkT , 
Pi - \ (5) 
Pgj , E B < QkT , 



where the parameter C, ~ 30 depends on magnetic field strength and temper- 
ature. 

The binding energy E B appearing in this formula is the energy required to 
release an ion from the surface of the condensed matter. This is equivalent to 
the energy required to release a neutral atom from the surface and ionize it, 
minus the e nergy gained by returning the electrons to the surface. Thus (see 
Tsong|ll990h . 



E B = Q s + Y,h-nct ) , (6) 



i=i 



where Q s > is the cohesive energy, Ii > is the iih ionization energy of 
the atom (so that Ya=i h the energy required to remove n electrons from the 
atom), n > is the charge of the ion, and > is the electron work function. 
All of these quantities were calculated in ML06b (see Figs. [T]and[2]). 



4-2 Conditions for Forming a Vacuum Gap 



If the ions or electrons that are required to fill the magnetosphere can do so 
freely, no gap will form. From Eqs. (j3J and ([SD, we can see that no polar gap 
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will form if 

<f> & C e kT ~ 2.6T 6 keV (7) 

for a negative polar magnetosphere (ft ■ B p > 0), and 

E B < dkT ~ 2.6T 6 keV (8) 

for a positive polar magnetosphere (fi • B p < 0). Note that the two conditions 
are related by almost the same factor of kT. This is coincidental, as C e and 
Ci depend on completely different constants and functions of B and T. 

For neutron stars in general, the electron work function is much less than 
30kT (see Fig. [2]), so electrons are easily removed from the condensed surface. 
No gap forms for a negative polar magnetosphere under neutron star surface 
conditions. The ion binding energy E B , on the other hand, is generally larger 
than 30kT. Ions can tightly bind to the condensed surface and a polar gap 
can form under certain neutron star surface conditions. Figure H] shows the 
critical temperature (determined by E B = CikT) below which a vacuum gap 
can form for the Fe, C, and He surfaces. 



5 Discussion 



Our calculations show that there are a range of neutron star magnetic field 
strengths and surface temperatures where the condensed surface will have an 
important effect on radiation from these stars. For example, if the surface com- 
position is Fe, then strong-field neutron stars (B > 10 13 G) with moderate 
(T <, 10 6 K) surface temperatures should have atmospheres/vapors that are 
effectively transparent to thermal radiation, so that the emission becomes that 
from a bare condensed surface. This may explain the nearly blackbody-like ra- 
diation sp ectrum observed from the nearby isolated neutron star RX J18 56.5- 
3754 (e.g. iBurwitz et aDl2003l : Ivan Adelsberg et alJlioOol ; IHo et alJlioofl ). 



Our work also indicates that for magnetar-like surface magnetic field strengths, 
a vacuum acceleration gap may form in the polar cap region of the neutron 
star, provided that the c ap temperature is not too h igh. As mentioned before, 



recent observations (e.g. jKaspi fe McLaughlinll2005l ) show that magnetars and 



radio pulsars overlap somewhat in their ranges of magnetic field strengths but 
most magnetars do not have pulsed radio emission. One possibility is that 
the magnetar also has radio emission, but the pulse is beamed away from 
us. On the other hand, magnetars typically have surf ace temperatures 3-5 
times larger than those of high-B radio pulsars (see, e.g. jKaspi fe Gavriilll2004 
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Fig. 4. The condition for the formation of a vacuum polar gap above He, C, and Fe 
neutron star surfaces, when the magnetosphere is positive over the poles (f2-B p < 0). 



and lKaspi &: McLaughlin! I2005I ). We therefore speculate that a key difference 
between magnetars and high-B radio pulsars is their difference in the surface 
temperature. If the surface binding energy is such that at the low surface 
temperature of a radio pulsar charged particles are held tightly while at the 
significantly higher surface temperature of a magnetar these particles flow 
freely out into the magnetosphere, under the polar gap model we will have the 
situation where the pulsar shows pulsed radio emission and the magnetar does 
not (see Fig. HJ). Note that the field strengths are inferred from the neutron star 
spin parameters, so it could be that these inferred strengths are incorrect. This 
would require a systematic error in the inferred field strengths of magnetars, 
due perhaps to the non-dipole nature of these superstrong magnetic fields. 



The r ecent detection of the radio emission from a transient AXP (ICamilo et al 



20061 ) is of great interest. The emission may have been triggered by a rear- 
rangement of the surface magnetic field, which made pair cascades possible. 
We note that the occurrence of pair cascades depends strongly on the field 
line curvature in the original vacuum gap model. 
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